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We review observations and gravitational lensing theory related to the magnification of back- 
ground QSOs by intervening overdensities, and the induced cross-correlation between sources and 
foreground galaxies. We pay special attention to simulations, and present some preliminary results 
from high resolution cluster simulations, which show the role of the halo core and substructure 
on non-linear magnification. For massive clusters, deviations from the weak gravitational lensing 
regime are significant on arcmin scales and bellow. The accumulated knowledge in the field already 
shows that gravitational lensing magnification is an important astrophysical and cosmological tool. 



OBSERVATIONS 

If we take two populations of astronomical objects sep- 
arated by a sufficiently large distance (for example galax- 
ies at low redshift and QSOs at large redshift) we expect 
that there will be no physical connection between them, 
and therefore the positions in the sky of the members of 
these two populations would be uncorrelated. However 
this is not the case in reality. 

In fact, several groups have been measuring the cross- 
correlation between the angular positions of objects at 
high redshits with objects at low redshifts. Bellow follows 
an incomplete, but representative, list. Some of the data 
sets are shown in Figures ^ and El (s refers to the double 
log slope of the cumulative number count as a function 
of flux for background luminous sources) . 

-1979, Seldner & Peebles Q found an excess of galaxies 
within 15 arcmin of 382 QSOs. 

-1988, Boyle, Fong & Shanks %^ (CS99 in the legend 
of Figure^ found an anticorrelation between faint high- 
redshift QSOs (s = 0.78) and low-redshift galaxies from 
machine measurements of photographic plates. 

-1997, Bemtez & Martmcz-Gonzales (BMG97 in 
Figure ^ found a positive cross- correlation between 144 
radio-loud PKS QSOs (s = 3.5) and COSMOS/UKST 
galaxies, and no correlation when using 167 optically se- 
lected LBQS QSOs (s = 2.5). 

-1988. WiUiams & Irwin % (WI98 in FigureCJ found a 
strong cross-correlation between optically selected QSOs 
(s = 2.75) from the LBQS Catalog and APM galaxies. 

-2001, Benitez, Sanz & Martmez-Gonzales % (BSMOO 
in Figure ^ found positive cross-correlations between 
radio-loud quasars from 1-Jy (s — 1.93) and Half-Jansky 
(s = 1.42) samples and COSMOS/UKST galaxies. 

-2003, Gaztanaga found a strong positive cross- 
correlation between QSOs (s = 1.88) and galaxies from 
the SDSS early data release. 

-2003 and 2005, Myers et al. H |i| found strong 
negative cross-correlations (anticorrelations) between ^ 
22, 000 faint 2dF QSOs (s = 0.725) and - 300, 000 galax- 
ies and galaxy groups from APM and SDSS early data 




FIG. 1: Compilation of some observational determinations of 
QSO-galaxy cross-correlation. The solid lines show theoret- 
ical predictions from analytical calculations assuming weak 
gravitational lensing, ACDM {Qm = 0.3, = 0.7, f2f, = 
0.019/^^, n = 1, h — 0.7, erg — 1.0), and foreground lens pop- 
ulations of power spectrum equal to the APM Galaxy survey 
(lower curve) and Abell-ACO Galaxy Cluster Survey (upper 
curve) . 

release. 

-2005, Scranton et al. 10] found cross-correlations be- 
tween ^ 200, 000 QSOs and ~ 13 million galaxies from 
the SDSS ranging from positive to negative signal, de- 
pending of the magnitude limits of the QSO population 
subsample (s = {1.95, 1.41, 1.07, 0.76, 0.50}). 

THEORY: ANALYTICAL 

The most successful hypothesis to explain the observed 
cross-correlations and anti-correlations is gravitational 
lensing. It generates two competing effects that can ex- 
plain both positive and negative cross-correlations be- 
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tween objects of two redshift distinct populations, in 
what is caUed magnification bias. 

The presence of a gravitational lens magnifies sources 
behind it, bringing to view sources that would be too 
faint to be detected in a magnitude-limited survey. This 
effect works to produce a positive cross-correlation be- 
tween objects physically associated to the foreground 
lenses and background objects that are magnified. On 
the other hand, the lens also enlarges the solid angle be- 
hind it, therefore the source density behind the lens is 
diluted, what works to produce a negative background- 
foreground cross-correlation. 

The factor that determines which of the two competing 
effects (magnification or dilution) is preponderant is the 
slope of the magnitude number count of the sources. If 
this slope is steep then many faint sources are brought 
to view, but if the slope is low (flatter) then few extra 
sources are brought to view and the dilution effects wins. 

To compute the cross-correlation between two popu- 
lations we calculate the correlation between the density 
contrast of the two groups 
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where Ug and Ug are the background and foreground pop- 
ulations (e.g. QSOs and galaxies or galaxy groups) den- 
sities. A bar over a quantity indicates its mean value, 
and (...) represents the average over cp and the direction 
of 6 (but not its modulus). This definition is equivalent 
to the cross-correlation estimator [llj 
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where DDiO) is the observed number of background- 
foreground pairs, and DR{9) is the expected number of 
random pairs. 

The ratio DD{6) / DR{0) is the enhancement factor due 
to the magnification bias, and under the assumption that 
the cumulative number counts by flux is of the form N{> 
S) oc S'"'* we have that 
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which cleary indicates that in an overdense region (/i > 
1): s > 1 leads to positive cross-correlation, and s < I 
leads to negative cross-correlation (anticorrelation) . 

The magnification can be written in terms of the grav- 
itational lensing convergence k and shear 7, 
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The magnification can be calculated in the cosmological 
context, assuming the Born approximation, as a weighted 
integration of the matter density field 

k(0)- W{y)5{e,y)dy , (5) 



where 5 is the density contrast, y is a comoving distance, 
yoo is the comoving distance to the horizon, and W{y) is 
a lensing weighting function 
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Gg is the source distribution, a is the scale factor, and 
fx is the curvature-dependent radial distance. 

The shear can be obtained from a convolution of the 
convergence j and therefore the knowledge of the mass 
distribution between the observer and the source plane 
allows the computation of the desired gravitational lens- 
ing effects. 

If we assume weak gravitational lensing, k <C 1, some 
analytical calculations become much simpler. The mag- 
nification becomes = 1 -I- 2k and the QSO -galaxy 
cross-correlation ^ can be expressed as UJI 
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where y is the comoving distance, which here parame- 
terizes time (2/00 represents a redshift z — 00), and k 
is the wavenumber of the density contrast in a plane 
wave expansion; Jo is the zeroth-order Bessel function 
of first kind; and fxiy) is the curvature-dependent ra- 
dial distance {— y for a flat universe). Pgm(k,y) can 
be seen as the galaxy- mass cross-power spectrum jisf . 
and under some assumptions Q| may be expressed as 
Pgm{k,y) = y/Pg{k)Prn{k,y), where Pg{k) is the power 
spectrum for galaxies or galaxy groups and Pm(k,y) is 
the non-linear time evolved mass power spectrum. 

Expression (0 indicates that the background- 
foreground cross-correlation due to lensing is depen- 
dent on several quantities of cosmological relevance. 
Guimaraes et al. 0| explores these cosmological de- 
pendences and Figure El illustrates the sensitivity of the 
cross-correlation between a population of QSOs at z=l 
and galaxies at z=0.2. 

Higher order terms can be added to the Taylor expan- 
sion of the magnification ^ producing a better approx- 
imation |l6j . However the full accounting of non-linear 
magnification is not feasible by this analytical path. 



SIMULATIONS 

The analytical approach has some limitations at fully 
incorporating deviations from the weak lensing approx- 
imation and properly modeling lensing selection issues. 
Therefore, computer simulations can be a useful com- 
plement for a better understanding of the problem. We 
describe, in the sequence, two kinds of simulations: one 




FIG. 2: Cross-correlation due to weak gravitational lensing 
dependence on cosmological model (top plot) and matter den- 
sity (bottom plot). The internal plots are the results for the 
mass power spectrum normalized to the cluster abundance 
(main curves use COBE normalization). Top plot: Solid lines 
are for SCDM, ^Im = 1, h = 0.5 (cts ~ 1.1); dashed ones for 
ACDM, = 0.3, Qa = 0.7, h = 0.7 [as = 1.0); and dotted 
ones for OCDM, Qm = 0.3, Qa = 0, h = 0.7 (erg = 0.46). 
Bottom plot: Dependence on matter density in a flat universe 
with cosmological constant {ilm + = 1). Solid lines are for 
flm = 0.6 {as — 1.4); dashed ones for £7™ = 0.4 (erg = 1.2); 
dotted ones for flm ~ 0.3 (erg — 1.0); and dot-dashed ones for 
fim = 0.2 (o-g = 0.72). Other parameters are fib = 0.019//i^, 
h = 0.7, n = 1 



aiming at very large regions and many lenses, collectively 
analyzed, but somewhat short on resolution on small 
scales, and another kind of simulation aiming at indi- 
vidual clusters and their substructure. 
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FIG. 3: Average magnification around mock galaxy groups. 
The number ranges in the legend are the number of galaxies 
for the sets of groups. Thick lines with filled symbols include 
departures from the weak lensing approximation, and thin 
lines with open symbols are the weak lensing approximation 
for the magnification calculation. Errors are the standard 
deviation of the mean. 



Galaxies and Galaxy Groups 

Both the mass and light (galaxy) distributions in the 
universe can be mocked from N-body simulations, and 
from those the cross-correlation between background and 
foreground objects due to gravitational lensing can be 
obtained 

The first step is to generate a representation of the 
mass distribution in a redshift cone from the observer 
(z = 0) to a source plane at high redshift. A galaxy mock 
catalog can be generated from the simulated density field 
and the adoption of a bias prescription for the galaxy 
population. This galaxy mock will have the galaxy den- 
sity and auto-correlation function desired, which can be 
set to mimic a chosen real galaxy survey. Also from the 
simulated density field the gravitational lensing for a cho- 
sen source plan can be calculated in the form of a lensing 
map (convergence, shear, or magnification) using the for- 
malism of the previous Section. 

Guimaraes et al. ^3 used the Hubble Volume Sim- 
ulation, a N-body simulation with 10^ particles of mass 
Mpart = 2.25 • lO^^/j-iMg in a periodic SOOO^/i-^Mpc^ 
box, initial fluctuations generated by CMBFAST, force 
resolution of 0.1 /i~^Mpc, and a "concordance model" pa- 
rameter set, = 0.3, = 0.7, T = 0.21, as = 0.90. 
The magnification map was calculated for a source plane 
at redshift 1, and the average magnification was mea- 
sured around chosen sets of lenses (galaxies or galaxy 
groups of varying membership identified in the projected 
sky mock) to determine the expected source-lens cross- 
correlation. 

Figure |3| shows the average magnification around 
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FIG. 4: Cross-correlation between QSO and mock galaxy 
groups. The number ranges in the legend are the number of 
galaxies for the sets of groups. Thick lines include departures 
from the weak lensing approximation, and thin lines are the 
weak lensing approximation for the magnification calculation. 
Curves for groups of membership 1 to 4 are shown individually 
(from bottom to top for low to high membership) . Errors are 
the standard deviation of the mean. 

galaxy groups, and Figure ^ shows the corresponding 
cross-correlation. Groups of larger membership trace 
denser regions, so have a higher magnification and 
stronger cross-correlation amplitude. 

Figure [S] compares the simulation results for galaxies 
and galaxy groups with observational data from Myers 
et al. d, Q ■ Simulation results for angles smaller than 1 
arcmin cannot be obtained due to limited simulation res- 
olution; however for angles from 1 to 100 arcmin the com- 
parison with data gives a large disagreement between the 
amplitudes of observed and simulated cross-correlations. 

This disagreement between some observed QSO-galaxy 
and QSO-group cross-correlations were already visible in 
Figure ^ and is source of controversy. 

High Resolution Galaxy Clusters 

One limitation of the simulations described in the pre- 
vious Section is the low resolution at small scales. To 
be able to probe small angular scales, and therefore re- 
gions near the cluster core or the halo substructure, it is 
necessary to use simulations of much higher resolution. 

Guimaraes et al. ^3 used high resolution cluster halo 
simulations carried out by the Virgo Consortium [T^ 
to study the gravitational lensing magnification due to 
galaxy cluster halos. These halos were generated by high 
mass resolution resimulations of massive halos selected 
from a large {A79h~^ Mpc)'^ N-body simulation. 

The magnification maps generated by the simulated 
cluster halos were calculated using equation |0J, assum- 
ing a source plan at redshift 1 and the cluster halo at 
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FIG. 5: QSO-galaxy and QSO-group cross-correlation. Ob- 
servational data is from Myers et al. 0,0 with field-to-field 
errors. Simulation uses groups with 9 or more galaxies and 
estimated error with same source density as observed data. 
Some observational points fall below the shown logarithmic 
scale. 

redshift 0.15. To evaluate the role of substructure it was 
also calculated the magnification of the homogenized halo 
in concentric rings, so the density profile is maintained, 
but the substructure is washed away. The weak gravita- 
tional lensing approximation was also used to calculate 
the magnification, so the departure from this regime can 
be quantified in the case of massive clusters. 

Figure ini shows for three cluster halos seen by three or- 
thogonal directions each, the average magnifications de- 
scribed above as a function of the angular distance to 
the halo center. Five other simulated clusters examined 
show similar curves (results not presented) . 

Departures from the weak lensing regime become im- 
portant at angles of few arcmin. At these same scales the 
contribution of substructure to the magnification can also 
be significant in some cases. 



DISCUSSION 

We reviewed some of the accumulated history of QSO- 
galaxy cross-correlation observations, and the gravita- 
tional lensing theory associated with it. 

The observations of cross-correlation in the sky be- 
tween populations of objects that are very apart in depth 
are old, numerous, and are becoming precise. The gravi- 
tational lensing explanation for the cross-correlation phe- 
nomenon is more recent and carries with it the possibil- 
ity of using these kind of measurements as a tool for 
cosmology and astrophysics. However, even though the 
magnification bias hypothesis is in qualitative agreement 
with observations in general, several measurements of 
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FIG. 6: Magnification curves for tliree simuiated ciusters viewed from tliree ortliogonal directions. Solid curves take in account 
the non-linear magnification due to cluster substructure; dashed curves are for the homogenized mass inside concentric rings; 
dotted curves use weak lensing approximation. 



the QSO-galaxy cross-correlation have a higher ampli- 
tude than what is predicted. 

The most recent and largest measurements of QSO- 
galaxy cross-correlation carried out using the approx- 
imately 200,000 quasar and 13 million galaxies of the 
Sloan Digital Sky Survey (SDSS) [13 are in agreement 
with theoretical predictions based on a "concordance 
model" and weak gravitational lensing (non-linear mag- 



nification was not taken into account by [lOj , but is con- 
sidered to be necessary by the authors for a more accurate 
modeling) . 

On the observational side, further work with the exist- 
ing data may clarify the reason for disagreements among 
different measurements of QSO-galaxy cross-correlation. 
If systematic errors are to blame for the huge amplitudes 
measured by various groups, then it is fundamental to 
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identify and characterize them. 

On the theoretical side, both analytical work and the 
use of simulations are helping to provide a realistic de- 
scription of the phenomenon. This theoretical framework 
in conjunction with observational data may be useful in 
determining quantities of astrophysical and cosmological 
interest, for example and most promisingly the average 
mass of lens populations and the galaxy-mass power spec- 
trum. 
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